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SUMMARY 


This  report  is  concerned  with  developing  mathematical  models  of  SSBN/ASW 
scenarios  which  can  be  used  to  compare  the  survivability  and  effectiveness  of 
present  and  future  strategic  force  mixes  against  a  variety  of  enemy  threats. 
The  work  was  supported  by  the  Office  of  the  Chief  of  Naval  Operations  through 
OP-604  under  the  direction  of  Mr.  Robert  Piacesi. 
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1.  INTRODUCTION 


The  policy  of  deterrence,  adopted  by  the  United  States  as  a  preventative  to 
the  outbreak  of  a  major  nuclear  war,  is  currently  embodied  by  a  mix  of  military 
forces:  long-range  bombers,  land-based  intercontinental  ballistic  missiles,  and 
submarine  launched  ballistic  missiles  (SLBM).  The  latter  member  of  the  mix  is  of 
particular  interest  because  it  is  generally  acknowledged  to  be  the  most  survivable 
against  an  all-out  attack.  Although  nuclear  submarines  are  virtually  undetectable 
by  present  technology,  a  major  problem  of  military  planning  is  estimating  submarine 
survivability  in  the  future,  and  the  consequent  effectiveness  of  the  surviving 
weapons  against  an  enemy  target  base.  The  purpose  of  this  report  is  to  develop 
some  mathematical  models  which  can  be  used  to  compare  both  survivability  and 
effectiveness  of  present  and  future  potential  SLBM  force  mixes  against  a  variety 
of  ASW  threats. 

The  principal  mathematical  tool  used  is  game  theory.  In  this  investigation, 
a  game  consists  of  exactly  two  moves,  one  by  each  of  two  opposing  players  labeled 
x  and  y.  The  x-player  moves  first,  followed  by  the  y-player  who  selects  his 
strategy  with  full  knowledge  of  his  opponent's  play.  Thus,  the  problem  of  the 
x-player  is  to  find  a  strategy  which  yields  minimal  advantage  to  the  y-player  when 
y  plays  his  best  response.  This  is  a  conservative  approach,  and  provides  a  "worst- 
case"  analysis  from  the  viewpoint  of  the  x-player. 

2.  SURVIVABILITY  MODEL 

We  begin  with  a  derivation  of  the  probability  distribution  used  to  describe 
the  survivability  of  a  single  submarine  which  comes  under  attack  at  time  t  *  0  by 
y  enemy  units  (ASW),  generally  taken  to  be  aircraft.  The  area  of  the  patrol  zone 
is  A.  The  aircraft  are  assumed  to  have  no  prior  knowledge  concerning  the  location 
of  the  submarine  within  the  zone,  and  therefore  search  for  it  randomly.  If 
detected,  the  submarine  is  destroyed.  The  motion  of  the  submarine  is  approximately 
stationary  compared  to  that  of  the  aircraft  over  a  small  instant  of  time  At,  so 
that  the  probability  of  detection  over  this  time  interval  is  SyAt/A,  where  S  is  the 
area  swept  out  by  one  aircraft  during  one  time  unit.  If  p(t)  denotes  the  probabil¬ 
ity  that  the  submarine  still  survives  at  lime  t,  then 


p(t+At) 


>(t+At)  - 
At 


-  ^r). 


-f  ■>“>• 


Letting  At-+0  yields  the  differential  equation 


The  initial  condition  p(0)  ■  1  determines  the  solution 


p(t)  - 
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We  are  now  ready  to  formulate  the  survivability  model  of  a  "surveillance- 
surge"  scenario.  The  Atlantic  and  Pacific  Oceans  are  each  partitioned  into  a 
number  of  patrol  zones.  Submarines  in  one  ocean  cannot  move  over  to  the  other, 
and  so  the  number  of  submarines  in  each  ocean  remains  constant.  Surveillance  does 
take  place  before  the  attack  begins,  reducing  the  amount  of  area  that  must  be 
searched  for  a  particular  submarine.  These  "areas-of-uncertainty"  have  been  calcu¬ 
lated  through  a  separate  simulation  procedure,  and  are  part  of  input  data.  The 
problem  is  to  determine  how  to  assign  submarines  to  patrol  zones  in  such  a  way  as 
to  maximize  the  number  of  weapons  surviving  an  attack  by  a  force  of  aircraft  which 
is  permitted  to  allocate  itself  optimally  amongst  specified  bases.  We  adopt  the 
notation: 

*  number  of  patrol  zones  in  ocean  k,  k  *  1,2 

-  number  of  submarines  assigned  to  ocean  k 

Y  *  total  number  of  available  aircraft 

-  area  of  uncertainty  for  a  submarine  in  patrol  zone  i 

T^  -  time  length  of  attack  in  zone  i 
w^  *  number  of  weapons  per  submarine  in  zone  i 
x^  ■  number  of  submarines  assigned  to  zone  i 
y^  *  number  of  aircraft  assigned  to  zone  i 

The  expected  number  of  surviving  weapons  is 


where  a^ 
aircraft. 


*1*2 


f(x,y) 


i«l 


WiXi* 


-Vi 


»  ST^/A^.  If  the  distribution  x  “  (x^jX^, • • • .x^  +N  )  is  known  to  the 
the  appropriate  response,  denoted  y(x),  is  determined  by 

f(x,y(x))  -  min  f(x,y). 

y 


Thus,  an  optimal  distribution  of  submarines,  x*,  must  satisfy 

f(x*,y(x*))  •  max  f(x,y(x))  -  max  min  f(x,y). 
x  x  y 

The  complete  statement  of  the  problem  is: 


max  min 
x  y 


12 

£ 

i-l 


WiV 


-Vi 


under  the  constraints: 
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(1) 

*1  +"‘+ “Bj  ‘  *1 

(2) 

\*i  +—+  X»1-m 

• 

(3) 

yi  +—+  ‘ 

(4) 

x±  >  0,  y,  ;  0 

•  X„ 
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Constraints  (1)  and  (2)  conserve  the  number  of  submarines  in  each  ocean,  and 
constraint  (3)  conserves  the  total  nuaber  of  aircraft;  constraint  (4)  eliminates 
meaningless  solutions  with  negative  assignments. 

3.  EFFECTIVENESS  MODELS 

Although  the  survivability  aodel  is  of  interest,  it  does  not  address  the 
principal  problem  of  estimating  the  daaage  capability  of  the  surviving  weapons  on 
an  enemy  target  base.  Because  its  sole  objective  is  survivability,  that  model 
might,  for  example,  predict  placing  most  submarines  far  from  the  target  base,  in 
effect  creating  an  arsenal  with  very  little  targeting  capability.  This  defect  is 
removed  by  incorporating  targeting  into  the  model  and  altering  the  payoff  function 
to  measure  the  expected  number  of  hard  targets  destroyed. 

A  base  of  targets  is  constructed  consisting  of  M  geographic  cells,  each  of 
which  contains  a  number,  v j ,  of  hard  targets.  The  probability  of  destroying  any 
target  within  cell  j  from  zone  i  is  assumed  to  be  equal  to  a  constant  value,  pjj. 
The  constraint  of  fratricide  requires  that  at  most  two  weapons  may  be  assigned  to 
any  one  target.  More  specifically,  varying  assumptions  on  second  weapon  assign¬ 
ments  give  rise  to  two  types  of  fratricide: 

(A)  the  second  laydown  weapon  assigned  to  a  target  may  be  launched  from  any 
zone  within  range; 

(B)  the  second  laydown  weapon  must  be  launched  from  the  same  zone  as  the 
first  weapon. 


The  fraction  of  surviving  weapons  in  zone  i  targeted  to  cell  j  in  the  first  laydown 

is  denoted  z^;  z^  is  given  the  analogous  meaning  for  the  second  laydown.  The 

order  of  the  player's  moves  is  an  extension  of  that  used  in  the  survivability  model. 
For  a  known  submarine  distribution,  the  aircraft  must  now  distribute  itself  in  such 
a  way  as  to  minimize  the  effect  of  any  possible  targeting  of  the  surviving  weapons. 
The  mathematical  statements  of  the  problems  corresponding  to  the  two  types  of 
fratricide  are: 
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M 

(A)  max  min  max  ^ 
x  y  z  j«l 


M 

(B)  max  min  max  £ 
x  y  z  j=l 

under  the  constraints: 


(1) 

X1  +**-+  \  *  X1 

(2) 

\+l  +*"+  "  X2 

(3) 

yl  +‘**+  yN1+N2  "  Y 

(4) 

_  v  _(1>  (2)  >  n 
ij*  yij‘  Zij  *  Zij  -  ° 

all  i,j 

(5) 

5  [•!!’  *  ■?] ’  ‘ 

i-1,2,"-,^^ 

(6) 

N.+N. 

v-2  _aiyi  (1)  < 

E  Vi®  zi.j  -  vj 

j“l»2, • • •  ,M 

Nl+N2 

E 

i«l 


_aiyi  (1) 
Vi6  PijZij 
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,  V 2  '  Vk 

+  L  wkv  p 

k»l 
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(7A) 

Nl+N2  -a,  y  Nl+N2  -ay 

V  Vk  (2)  <  V  aiyi  (1) 

£  Vk®  Zkj  -  E  Vi®  zij 

J*l»2, • • •  ,M 

(7B) 

zj2)  <  z<f  all  i,j 

Constraint  (5)  prevents  assigning  more  weapons  than  are  surviving  in  zone  1, 
constraint  (6)  limits  the  number  of  weapons  assigned  to  cell  j  in  the  first  laydown 
to  be  no  more  than  the  number  of  targets  within  that  cell,  and  constraints  (7A)  and 
(7B)  are  the  mathematical  descriptions  of  the  two  types  of  fratricide  corresponding 
to  the  payoff  functions  (A)  and  (B). 
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Related  models  occur  when  the  players  measure  their  respective  payoffs 
differently.  Of  particular  interest  is  the  case  in  which  the  x-player  maximizes 
the  number  of  hard  targets  destroyed,  and  the  y-player  minimizes  the  surviving 
number  of  weapons.  Mathematically,  this  problem  is  identical  to  the  one  described 
in  the  preceding  paragraph  except  that  the  payoff  function  of  the  y-player  is  now 
equal  to  the  payoff  used  in  the  survivability  model. 

4.  SOLUTIONS 


The  algorithms  discussed  in  this  section  provide  continuous,  rather  than 
discrete,  answers.  Because  the  discrete  problems  are  very  large  and  complex,  this 
has  become  a  generally  accepted  practice  in  dealing  with  these  types  of  models. 

The  continuous  answers  are  sometimes  given  probabilistic  interpretations. 

An  algorithm,  based  on  the  notion  of  hedging,  is  used  to  solve  the  survivabil¬ 
ity  problem.  We  visualize  a  series  of  plays  of  the  game  in  which  the  players 
respond  to  each  other  alternately.  Suppose  the  x-player  Initiates  play  with  the 
strategy,  xq.  The  y-player  responds  optimally  with  y^  determined  by 

fCXg.y^  -  min  f(xQ,y). 

y 

The  x-player  counters  against  y^  with  x^  determined  by 

f^.y^  ■  max  f (x,y^) , 
x 

and  the  y-player  likewise  responds  with  y«.  Now  the  x-player  hedges  against  yj  and 
y2,  that  is,  he  plays  xj  determined  by 

min  f(x,,y.  )  -  max  min  f(x,y.  ). 
k-1,2  K  x  k-1,2 

The  play  alternates  in  an  analogous  manner  with  xn  and  yn  satisfying 

min  f  (x^y^)  "  max  mip  f(x,yk) 
k-l,***,n  x  k-l,***,n 


and 

f(Vyn+l>  "  nin  f(Vy)’ 

y 

respectively.  If  x*  is  a  limit  point  of  {x  }  and  y*  is  the  optimal  ASW  response  to 
x*,  then  (x*,y*)  is  a  saddle  point  solution?  For  a  complete  theoretical  discussion 
of  hedging  and  its  application  to  the  survivability  model,  see  Winston  [1]. 


1.  Winston,  E.,  "Hedging  and  Maxmln,"  submitted  for  publication. 
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For  a  fixed  submarine  distribution,  the  hedging  algorithm  is  applicable  to  the 
inner  "min  max"  problem  of  the  effectiveness  models.  In  doing  so,  however, 

y  * 

modifications  are  made  in  determining  the  solutions  of  each  of  the  hedges.  The 
targeting  portion  of  the  algorithm  is  a  very  large  linear  programming  problem 
which,  if  solved  directly,  would  result  in  a  prohibitive  time  and  storage  require¬ 
ment  for  finding  a  solution  of  a  "min  max"  problem.  Therefore,  a  marginal  analysis 

y  z 

procedure,  based  on  an  ordering,  in  magnitude,  of  the  numbers  {p^.p^Cl-p  )},  Is 

used.  Moreover,  the  targeting  solutions  determined  by  this  method  automatically 
satisfy  all  the  constraints.  On  the  other  hand,  constraints  (6)  and  (7A)  are  non¬ 
linear  in  y  and  so  are  relaxed  when  finding  the  ASW  hedging  solutions.  Although 
not  rigorously  optimal,  this  approximation  is  very  fast  computationally  and  allows 
the  hedging  algorithm  to  converge  to  good  answers. 

Because  it  is  not  possible  to  determine  a  closed  form  for  F(x)  - 
min  max  f(x,y,z),  a  search  method,  based  on  a  scheme  proposed  by  Box  [2],  is  used 

y  z 

to  maximize  F(x).  F(x)  is  evaluated  over  a  random  sample  of  size  of  feasible 

points,  and  the  worst  point  is  identified.  The  centroid  of  the  remaining  points  is 
feasible  because  the  constraints  are  linear  in  x.  It  follows  that  every  point  with 
non-negative  components  which  lies  on  the  line  determined  by  the  worst  point  and 
the  centroid  is  also  feasible.  A  point  which  is  a  times  as  far  from  the  centroid 
as  the  reflection  of  the  worst  point  through  the  centroid  is  first  tested  for 
feasibility  and  then  compared  to  the  worst  point.  This  point  is  contracted  half¬ 
way  to  the  centroid  if  it  is  neither  feasible  nor  an  improvement  over  the  worst 
point.  If  an  admissible  point  is  found  in  fewer  than  c  contractions,  it  replaces 
the  worst  point,  and  the  procedure  is  repeated;  otherwise  the  best  point  is 
retained  and  Nj+^-l  new  random  points  are  generated.  The  process  stops  when  the 
distance  between  a  centroid  and  replacement  point  is  sufficiently  small. 

5.  SOFT  TARGETING 

The  effectiveness  models  developed  in  Section  3  have  a  target  base  consisting 
of  hard  targets  only.  Naturally,  the  expected  number  of  hard  targets  destroyed  is 
reduced  when  other  types  of  targets  are  also  taken  into  account.  This  phenomena  is 
observed  when  a  soft  target  base,  consisting  of  conventional  military  targets  and 
economic  targets,  is  incorporated  into  the  model.  Because  the  amount  of  hard 
target  damage  remains  as  the  primary  objective,  specified  amounts  of  soft  target 
damage  are  inserted  as  constraints.  Marginal  analysis  is  performed  on  soft  target 
damage  functions  of  the  form 


d(x)  ■  ve 


2.  Box,  M.  J. ,  "A  New  Method  of  Constrained  Optimization  and  a  Comparison  to  other 
Methods,"  Computer  Journal,  8,  1965,  pp.  42-52. 
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where  v  is  the  target  value,  and  H  is  a  fitting  parameter.  The  .first  weapons  used 
in  soft  targeting  are  those  surviving  weapons  which  have  not  previously  been 
assigned  to  hard  targets.  When  this  stock  has  been  depleted,  weapons  are  removed 
from  hard  targets  in  the  reverse  order  of  their  assignment,  and  reassigned  to  the 
appropriate  soft  targets. 

6.  CONCLUSION 


In  this  report,  a  number  of  game-theoretic  models  of  various  "surveillance- 
surge"  scenarios  have  been  presented.  The  most  interesting  and  pertinent  models 
contain  the  following  key  features: 

(1)  maximization  of  hard  target  damage, 

(ii)  inclusion  of  a  soft  target  base, 

(ill)  consideration  of  fratricide  constraints. 

Computer  codes  for  all  models  have  been  developed  and  are  being  exercised. 
The  results  of  a  series  of  numerical  studies  will  appear  in  the  classified 
literature. 
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